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Abstract. We estimate the distance modulus to long gamma-ray bursts (LGRBs) using 
the Type I Fundamental Plane, a correlation between the spectral peak energy E p , the peak 
luminosity L p , and the luminosity time Tl (= E- lso /L p where E- lso is isotropic energy) for small 
Absolute Deviation from Constant Luminosity(ADCL). The Type I Fundamental Plane of 
LGRBs is calibrated using 8 LGRBs with redshift z < 1.4. To avoid any assumption on the 
cosmological model, we use the distance modulus of 557 Type la supernovae (SNela) from the 
Union 2 sample. This calibrated Type I Fundamental Plane is used to measure the distance 
moduli to 9 high-redshift LGRBs with the mean error = 0.31, which is comparable with 
that of SNe la = 0.26 where [i stands for the distance modulus. The Type I Fundamental 
Plane is so tight that our distance moduli have very small uncertainties. From those distance 
moduli, we obtained the constraint S^m = 0.22 ± 0.04 for flat ACDM universe. Adding 
9 LGRBs distance moduli (z > 1.4) to 557 SNela distance moduli {z < 1.4) significantly 
improves the constraint for non-flat ACDM universe from (Om,^a) = (0.29±0.10, 0.76 ±0.13 
) for SNela only to (fl M , O\)=(0.23 ± 0.06, 0.68 ± 0.08) for SNela and 9 LGRBs. 
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1 Introduction 

Since the discoveries of possible distance indicators of long gamma-ray bursts (LGRBs) [1— 
7], many attempts to measure the distance up to the high redshift universe (z > 1.7) which 
cannot be probed by Type la supernovae have been done to test cosmological models [8-20]. 
These attempts are very important since no other distance indicators have been known in 
such high redshift universe. 

However the large statistical and systematic errors of the distance indicators restricted 
the precision of distance measurements by LGRBs. To estimate distances of LGRBs more 
precisely, there are two directions : to discover more precise distance indicator [6, 10, 21] 
and/or to take weighted average of distances derived by many distance indicators [11, 13, 16]. 
These studies, however, are on the verge of many difficulties, such as missing or multiple jet 
breaks in Swift era [22, 23] or increasing scatters of distance indicators [24-26]. To make 
LGRBs to be more precise and reliable standard candles, we must know the origins of the 
dispersion of distance indicators to control them carefully. 

In our previous studies, we have investigated origins of dispersion of the spectral - bright- 
ness correlations and found possible origins of systematic errors [27, 28]: (i) spectral peak 
energies (E p ) estimated by fitting the spectra with cut-off power-law model with three free 
parameters tend to be overestimated compared with those obtained by using the Band model 
with four free parameters [29], (ii) peak luminosities (L p ) should be estimated for a fixed 
time scale in the individual rest frame rather than the observer frame, (iii) contamination of 
short GRBs or other unknown populations. Removing the systematic errors by (i) and (iii), 
and correcting (ii), we found that the Fundamental Plane of LGRBs, that is, a correlation 
between the spectral peak energy E p , the rest frame 2.752-second peak luminosity L p , and 
the luminosity time Tl (= E- lso /L p where E{ so is the isotropic energy), is much tighter than 
the original E p — Tl — L p correlation by R. Tsutsui et al. [17] [See also 30]. 

Recently, we found that dividing LGRBs into two groups according to the absolute 
deviation from their constant luminosity (ADCL) significantly improves the Fundamental 
Plane of LGRBs [31]. ADCL is of critical importance in this work so let us summarize the 
main results below. ADCL is a quantity which characterizes the overall shape of light curve 
of a GRB and defined as follows: 




(1.1) 
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where Cf orm is normalized cumulative counts and tf orm is the normalized time by Tgg in 
which the cumulative counts increase from 0.01 to 0.99. Most of small- ADCL(< 0.17) events 
consist of many pulses with almost the same brightness, while most of large- ADCL{> 0.17) 
events consist of a bright main pulse and the weak long tail like FRED (Fast Rise and 
Exponential Decay). In Figure. 1, we show the cumulative light curves of LGRBs from 
Tsutsui, et al. [31]. The time and the cumulative counts are normalized by Tgs and the total 
counts, respectively. The red solid lines indicate small- ADCL events while blue solid lines 
do large- ADCL events. The black dashed line indicates the line of the virtual source with a 
constant luminosity. We clearly see that red lines with small ADCL (< 0.17) and blue solid 
lines with large ADCL(> 0.17) form different groups in this figure so that our classification 
scheme for subclasses of LGRB would be reasonable. 

In Tsutsui, et al. [31], we found that the best-fit function for small ADCL(< 0.17) 
events is given by 

/ p \ 1.84±0.03 
L 1() 52.53±0.01 ergs -l ( 



10 2 - 71 keV, 

/ rp \ 0.29±0.08 

with xt/d-o.f = 10.93/14 and o"i n t = 0, while that for large ADCL(> 0.17) events is given by 

/ rp \ 1.82±0.093 
-,52.98±0.080 / -^P 



Lp ~ 10 6rgS Vl0 2 - 71 keV 

/ rj-\ \ 0.85±0.26 

with xl/d-o.f = 7.58/8 and o"int = 0. Since both best-fit functions separately form the plane 
in (log E p , log Tl, log L p ) space, we call the best fit function for small-A-DCL (large- ADCL) 
as Type I (Type II) Fundamental Plane, respectively. Figure. 2 shows the Type I (left) and 
Type II (right) Fundamental Planes for given cosmological parameters from Tsutsui, et al. 
[31]. In both figures, Small- ADCL events are marked with red filled circles and large- ADCL 
events are with blue squares. Orange and light blue squares indicate outliers from the Type 
I and Type II Fundamental Planes, respectively, (see [31] for the definition of outliers. ) 
The left (right) of Figure. 2 is the plane perpendicular to the Type I (Type II) Fundamental 
Plane so that the Type I (Type II) Fundamental Plane is expressed by the solid line in each 
figure. In the left of Figure. 2, blue triangles which belong to the large ADCL events are 
above the solid line (=the Type I Fundamental Plane) while in the right of Figure. 2, red 
filled circles which belong to the small ADCL events are below the solid line (=the Type II 
Fundamental Plane) so that the existence of two Fundamental Planes can be recognized. The 
existence of such two subclasses is very similar to that of two Period - Luminosity relations 
of Cepheid variables [32]. 

In this paper using 21 small-AL>CL LGRBs from the 36 gold sample of Tsutsui et al. 
[31], we try to measure distance modulus up to z = 3.57 and determine the cosmological 
parameters. To avoid any assumption on cosmological models, we use the distance moduli of 
557 Type la supernovae (SNela) from the Union 2 sample to calibrate the Type I Fundamental 
Plane of LGRBs by the 10 low-z events in 0.168 < z < 1.30. This is a big difference from 
our previous paper [31] where cosmological parameters were assumed to derive the relation. 
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Figure 1. The cumulative light curves of LGRBs from Tsutsui, et al. [31]. The red solid lines indicate 
small- ADCL events while blue solid lines do large- ADCL events. The black dashed line indicates 
the line of the virtual source with the constant luminosity. We clearly see that red lines with small 
ADCL (< 0.17) and blue solid lines with large ADCL(> 0.17) form different groups. 



Applying the Type I Fundamental Plane calibrated at low redshifts without cosmological 
parameters to high-z GRBs, we obtain the relation between the redshift and the distance 
modulus \jl of 11 high-z events in 1.54 < z < 3.57. From these distance moduli, we will 
obtain the constraint on J1m = 0.22 ± 0.04 for flat A CDM universe. In this paper we do 
not use the Type II Fundamental Plane as a distance indicator to constrain the cosmological 
parameters due to the lack of enough number of large ADCL events. Throughout the paper, 
we fix the Hubble parameter as Hq = 70 km s _1 Mpc . 

2 Calibration 

2.1 local regression 

First, we calibrate the Type I Fundamental Plane of LGRBs using low-z events whose distance 
moduli can be obtained by those of Type la supernovae. Actually, to estimate the distance 
moduli of low-z GRBs, we use a local regression method first developed by W. S. Cleveland. 
[33] and W. S. Cleveland and S. J. Devlin. [34], and first applied for GRB cosmology by V. F. 
Cardone et al. [16]. This method does not need a global shape function of redshift-distance 
modulus relation and is one of interpolation procedures which use only the data near points 
and fit them with a low-degree polynomial. Although the number of data and the degree of 
polynomial are arbitrary variables, it does not depend on a shape of a global function. 

Applying this method, we obtain the distance moduli of 10 low-redshift LGRBs with 
mean error = 0.065 where /i stands for the distance modulus. 
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Figure 2. The E p -Tl-L p diagram both for the Type I (left) and Type II (right) Fundamental Planes 
for given cosmological parameters from Tsutsui, et al. [31]. Small- ADCL events are marked with red 
filled circles and large- ADCL events with blue squares. Except for outliers ( orange triangles), most 
of small- ADCL events are on the Type I Fundamental Plane (=thc solid line in the left figure) within 
their error bars without any systematic errors, while most of large- ADCL events are on the Type II 
Fundamental Plane (=the solid line in the right figure) except for outliers (light blue triangles), (sec 
[31] for the definition of outliers. ) In the left figure, blue squares which belong to large ADCL events 
are above the solid line (=the Type I Fundamental Plane) while in the right figure, red filled circles 
which belong to small ADCL events are below the solid line (=the Type II Fundamental Plane) so 
that the existence of two Fundamental Planes can be recognized. The existence of such two subclasses 
is very similar to that of two Period - Luminosity relations of Cepheid variables. 



2.2 calibration 

Having estimated the distance moduli of 10 low-redshift GRBs in a model independent way, 
we can calibrate the Type I Fundamental Plane of GRBs. Let us assume a linear correlation 
between E p , Tl and L p in logarithmic scale as, 

logT p (T p ,T L ) = A + Blog {E p /E p ) + Clog (T L /T L ) , (2.1) 

where A, B and C are the parameters of the model and E p and Tl are logarithmic mean 
of E p and L p with 10 2 ' 43 keV and 10 a97 sec, respectively. From equation (2.1), we obtain 
H(E P ,T L ,F P ) given by, 

(i(E p ,T L ,F p ) = 2.51ogT p (£ p ,T L ) - 2.5 log(4^F p ) 

-51og(3.0857 x 10 19 ), (2.2) 

where F p is the observed energy flux. We define a chi square function, 

X \A, B,C)ocj: ^;^ E f L 2 ,Fp))2 ' (2-3) 

1=1 °rtocal>« ""VgRBi* 
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and 



.2 

MGRB,* 



T [(i + 2cK gFp4 

_L d2^2 i ,^2^2 | 2 

+ B a \ogE v;i + C log Tl,i + ^ir 



int 



(2.4) 



where o"i n t is the systematic error determined by iteration. 

Before we minimize the chi square function, we use the robust statistic and outlier de- 
tection method which have been developed in Tsutsui et al. [30]. As a result, we detected 
two outliers (091003, 080319B) , and removed them from following chi square analysis. Per- 
forming chi square regression on the remaining samples, we obtained the best fit parameters 
and l-o" uncertainties summarized in table. 1. The parameters in table. 1 and outliers of 
the relation are consistent with those in Tsutsui et al. [31]. Figure. 3 shows the Ep-T^-Lp 
diagram for small- ADCL events in z < 1.4. Two outliers are marked with green triangles 
and the solid line indicates the best fit function of Eq.(2.1) with values of A, B and C in 
table 1. 



Figure 3. The E p -Tl-L p diagram for small- ADCL events (red filled circles) in z < 1.4 without any 
assumption of cosmological parameters. The solid line is the best fit function with xt = 0-83 and 
0mt = 0. The green triangles are outliers, (sec [31] for the definition of outliers.) 



Table 1. The best fit values (1-cr uncertainties) of the model parameters of the Type I Fundamental 




10" 



cr 1 10° io 1 



(E p /l0 2 - 43 [keV]) 19 (T L /l0 097 [sec])' 



,0.25 



Plane (Eq. (2.1)). 




52.06(0.02) 1.90(0.05) 0.25(0.13) 0.83 
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Table 2. The redshifts, distance modulus of LGRBs estimated by local regression method and the 
ones estimated from the Type I Fundamental Plane. The first through the seventh column express 
name of the GRB, the redshift, the distance modulus determined by Typela supernovae, the error of 
Mlocaij the distance modulus determined by the best fit function of Type I Fundamental Plane (Eq. 
(2.1)), the error of /z, the value of ADCL, respectively. Four outliers are separated in the last four 
rows. They are the same as the outliers for small- ADCL GRBs in Tsutsui et al. [31] 



GRB 


z 


/^local 


r^iocai 


f 1 




ADCL 


030329 


0.168 


39.48 


0.04 


39.61 


0.17 


0.13 


090618 


0.54 


42.43 


0.04 


42.33 


0.19 


0.10 


050525 


0.606 


42.83 


0.07 


42.72 


0.13 


0.13 


970228 


0.695 


43.15 


0.05 


42.7 


0.72 


0.14 


990705 


0.843 


43.75 


0.06 


43.86 


0.22 


0.10 


091208B 


1.063 


44.3 


0.09 


44.77 


0.77 


0.07 


061007 


1.261 


44.73 


0.06 


44.78 


0.16 


0.14 


990506 


1.3 


44.81 


0.07 


44.82 


0.26 


0.08 


070125 


1.547 


- 


- 


45.39 


0.25 


0.10 


990123 


1.6 






45.5 


0.16 


0.09 


990510 


1.619 






45.46 


0.34 


0.14 


090902B 


1.822 






45.91 


0.15 


0.09 


081121 


2.512 






46.9 


0.24 


0.07 


080721 


2.602 






46.6 


0.28 


0.16 


050401 


2.9 






46.5 


0.42 


0.15 


971214 


3.418 






47.5 


0.73 


0.12 


090323 


3.57 






48.05 


0.3 


0.12 


091003 


0.897 


43.82 


0.09 


45.97 


0.28 


0.16 


080319B 


0.937 


43.82 


0.08 


45.36 


0.17 


0.14 


110213A 


1.46 






42.9 


0.38 


0.09 


110422A 


1.77 






43.51 


0.16 


0.10 



3 Constraint 



Applying equation (2.2) to high redshift LGRBs, we can obtain the extended Hubble diagram 
with the unprecedented precision. We define cr Mi as below, 

n 25 



a 



a\ + (a B log E p /E p y + (a c log T L /T L ) 



+ (Ba logEp ) 2 + (Ca logTL ) 2 + (1 + 2C> 1 2 ogFp 



(3-1) 



hit- 



Here we neglect the effect of the gravitational lensing because M. Oguri and K. Takahashi. 
[35] discussed the gravitational lensing of GRBs and found that the biases are small. 

Then we derive constraints on cosmological parameters. In the A-CDM model, the 
luminosity distance [di) is given as a function of the density parameters, VLm and Oa- 

The chi square function is defined by, 



N 

x 2 (n M ,n A ) = ]T 
i=i 



Hi - // h (zi,f2 M ,^A,) 



-, 2 



(3.2) 
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where [i th (zi, 0*m, ^a) = 51og(d^ 1 /Mpc) + 25. We used the robust regression and outlier 
detection method again. As a result, we detect two outliers (110213A, 110422A), and removed 
it from following chi square analysis. In table 2, we summarize the redshifts, distance modulus 
of LGRBs estimated by local regression method and the ones estimated from the Type I 
Fundamental Plane with their 1-a uncertainties. 




Figure 4. Extended Hubble diagram from the Type I Fundamental Plane. Gray squares: the 557 
SNe la from R. Amanullah et al.[36]. Blue open circlues: LGRBs at z < 1.4 estimated from the local 
regression of SNe la data, red filled circles: LGRBs estimated from the Type I Fundamental Plane. 
Green triangles: outliers. The solid line indicates a theoretical model calculated with (Om, fiA)=(0.30, 
0.70). Inset : A residual Hubble diagram from the theoretical model. 



Figure 4 shows an extended Hubble diagram up to z = 3.57. The gray squares in- 
dicate the SNe la from R. Amanullah et al. [36], the blue open circles indicate the low 
redshift LGRBs at z < 1.4 estimated from the local regression of SNe la data. Red 
filled circles and green triangles indicate LGRBs estimated from the Type I Fundamental 
Plane, and outliers, respectively. The solid line indicates a theoretical model calculated with 
(r^Mj ^a)=(0.30,0.70). As Figure 4 shows, the uncertainties of distance modulus of LGRBs 
are as small as those of individual SNe la (gray squares). The mean error of distance modulus 
of LGRBs are a^ = 0.31, while those of SNe la are a^ = 0.26. 

Performing the chi square analysis on the remaining 9 LGRBs, we obtained the con- 
straint Qm = 0.22 ±0.04 for flat ACDM universe. In Figure 5, we show the likelihood contour 
for non-flat ACDM universe from 9 LGRBs (green), with the likelihood contours from 557 
SNe la (blue) and 9 LGRBs + 557 SNe la (red), and the best-fit values with 1-a errors are 
shown in table 3. Figure 4 shows that adding only 9 high redshift LGRBs (z > 1.4) to 557 
SNe la {z < 1.4) significantly improves the constraint on f2jy[ _ ^A plane. 
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0.2 0.4 0.6 0.8 1 

Figure 5. Constraints on (£Im, ^a) plane from 9 LGRBs with z > 1.4 (green), 557 SNe la (blue) and 
9 LGRBs + 557 SNe la (red). The contours correspond to 68.3%, 95.4%, 99.7% confidence regions. 
See also Table 3. 







n A 


xVd.o.f 


9 LGRBs (flat) 


0.22 ± 0.04 




5.09/7 


557 SNe la (non-flat) 


0.29 ±0.10 


0.76 ±0.13 


542.1/555 


Combined (non-flat) 


0.23 ± 0.06 


0.68 ±0.08 


548.3/564 



Table 3. Constraints on (Qm, ^a) for flat universe from 9 LGRBs and for non-flat universe from 557 
SNe la and 9 LGRBs + 557 SNe la with their reduced chi squares. All the quoted errors are at the 
1-cr confidence level. 



4 Discussion 

In this paper we extended the Hubble diagram with the Type I Fundamental Plane of LGRBs. 
Because of the tightness of the relation, our distance measurements are much more precise 
than those of previous works. Although the number of LGRBs are much smaller than that of 
the previous works [11-17], our constraint is much stronger. This indicates the importance of 
controlling and removing systematic errors. Although we used only the Type I Fundamental 
Plane to measure distance modulus, in principle, we may also use the Type II Fundamental 
Plane as an independent distance indicator. However, because there is still a small number of 
long tailed events, it is difficult to calibrate the relation with our outlier rejection technique. 

The constraints from 9 high redshift GRBs are consistent with those from other probes, 
e.g. Cosmic Microwave Background (CMB), Baryon Acoustic Oscillation (BAO)[37, 38]. 
We should stress here that this consistency is not trivial. Since the redshift range of GRB 
(1.4 < z < 3.5) is different from Typela supernovae (z < 1.7), BAO (z ~ 0.35) and CMB 
[z = z rec ), they might show different cosmological parameters if the dark energy is strongly 
time dependent. Therefore our results suggest that time dependence of the dark energy is 
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not so strong even if it exists. 

The purpose of this paper is to show the potential high ability of Type I Fundamental 
Plane of the LGRBs as a distance indicator. Combining our data with other cosmological 
probes is the next step. A different direction is to test more general equation of state than 
ACDM universe. We are now trying to combine our data with CMB and BAO as well as 
more general equation of state to constrain the property of the dark energy. 
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